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ABSTRACT 

We present the first systematic investigation into the metallicity gradients in galaxy close pairs. We 
determine the metallicity gradients for 8 galaxies in close pairs using H II region metallicities obtained 
with high signal-to-noise multi-slit observations with the Keck LRIS Spectrograph. We show that the 
metallicity gradients in close pairs are significantly shallower than gradients in isolated spiral galaxies 
such as the Milky Way, M83, and M101. These observations provide the first solid evidence that 
metallicity gradients in interacting galaxies are systematically different from metallicity gradients in 
isolated spiral galaxies. Our results suggest that there is a strong relationship between metallicity 
gradients and the gas dynamics in galaxy interactions and mergers. 

Subject headings: galaxies:starburst — galaxies: abundances — galaxies:fundamental parameters — 
galaxies:interactions 



1. INTRODUCTION 

Galaxy interactions and mergers are fundamental to 
the formation and evolution of galaxies. Current N- 
body/smoothed particle hydrodynamic (SPH) simula- 
tions predict that, as a merger progresses, galaxy disks 
become disrupted by tidal effects, causing large gas in- 
flows into the cent ral regions where kpc-scale starbursts 
may be trigge red dBarnes fc Hernauistl l!996|). or AGN 
may be fueled (Sori ngel et al.ll2005HHopkins et al.1 12006). 
At very late stages, starburst-or AGN- driven super- 
winds may drive gas back out of the central regions 
(|Naravanan et al.ll2008l) . 

Evidence for major gas inflows in i nteracting galax- 
ies h as been seen in ionized gas (iRampazzo et all 
20051). in neutral gas dHibbard fc van Gorkoml If 9961: 
Georgakakis et all 120001: l lono et aLl l2005t lEmonts et al.l 
20061 : ICuUen et al.l 120071) . and in kinematic studies us- 



ing absorption lines (|Rupke et al.l I2UU51; iMartinl l2uutl . 
Whether large-scale gas flows occur may depend on the 
mas s difference between the interacting galaxies (see 
e.g..lWoods et al.ll2006l : iHorellou fc KoribalskilfeOOTl ) . In 
iKewlev et al.l (|2006aD . we compared the nuclear gas- 
phase metallicities of nearby field galaxies with nearby 
galaxy pairs. We found that close pairs have systemati- 
cally lower metallicities than either field galaxies or more 
widely separate d pairs at the same luminosity (see also 
iLee et al.1 120041 ). Most isolated late-type spiral galaxies, 
and starburst barred spirals display strong metallicity 
gradients (see Henry & Worthey f999 for a review). We 
proposed a merger scenario in which galaxy interactions 
drive large gas flows towards the central regions, carry- 
ing less enriched gas from the outskirts of the galaxy into 
the central regions, disrupting metallicity gradients and 
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diluting central metallicities. 

Subsequent research into central metallicities of inter- 
acting galaxies confirms the relationship be tween galaxy 
intera ctions and central metallicities. Rupk eet al.l 
(120081 ) show that merging luminous infrared galaxies 
have lower central metallicities compared with local iso- 
lated galaxies of similar luminosity or mass, suggest- 
ing that large interaction-induced gas flows have diluted 
the centr al metallicity of lu mino us infrared galaxies. 
Simila rly, lEllison et ail (|2008| ) and iMichel-Dansac et al.1 
(2008) showed that galaxy pairs in the Sloan Digital 
Sky Survey have lower central metallicities for their 
stellar mass than isolated galaxies. Outliers from 
the standard mass-metallicity relation for star forming 
galaxies have been shown t o be morphologica lly dis- 
turbed jPesElesiran HH; EonsoeLan [2013) . Ad- 
ditional investigations suggest a relation ship between 
large - scale environment and metallicity (|Cooper et al.l 
2008 UEllison et al.l [20091) . 

iChien et al.l (|2007| ) investigated the stellar population 
ages and metallicities of 12 young star clusters in the 
merging galaxy pair NGC 4676. They found a relatively 
flat distribution of oxygen abundances along the north- 
ern tail, suggesting e fficie nt gas mixing within the tail. 
iTrancho et all (|2007aD and lBastian et all (|2009l) find sim- 
ilarly flat metallicity distributions in star clusters in the 
merging galaxies NGC 3256 and the Antennae, respec- 
tively. 

Theoretical simulations that include the effect of merg- 
ers on the metallicity distribution in galaxies support the 
picture of major gas inf lows induced by the t idal effects of 
galaxy interactions. In lRupke et all (|2010af ). we used N- 
body/SPH simulations to investigate the theoretical rela- 
tionship between metallicity gradients and galaxy major 
mergers. These simulations are based on close passage, 
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TABLE 1 

Sample Galaxies 



ID 


Common 


Coordinates 


M B 


Nuclear 


Hubble 


Separation 


Central 


R25 


Gradient 




Name 


(J2000) 




Class a 


Type b 


(kpc) c 


log(0/H)+12 a 


(kpc) 


Slope (R/R 25 ) d 


1 


UGC 12914 


00 01 38.3 +23 29 01 


-20.88 


HII 


S(r)cd? pec 


17.4 


9.08 


14.69 


-0.104 ±0.048 


2 


UGC 12915 


00 01 41.9 +23 29 45 


-20.09 




S? 


17.4 




11.11 




3 


UGC 312 


00 31 23.9 +08 28 01 


-20.53 


HII 


SB? 


21.1 


8.78 


11.67 


-0.133+0.079 


1 


UGC 813 


01 16 16.4 +46 44 25 


-20.44 


HII 


Sb 


16.5 


9.00 


11.16 


-0.233+0.045 


5 


UGC 816 


01 16 20.5 +46 44 53 


-20.78 


HII 


S? 


16.5 


8.91 


13.33 


-0.252+0.032 


6 


NGC 3994 


11 57 36.9 +32 16 39 


-19.85 


C 


SA(r)c pec 


24.7 


8.93 


6.18 


-0.139+0.045 


7 


NGC 3995 


11 57 44.1 +32 17 39 


-20.67 


HII 


SAm pec 


24.7 


8.75 


19.40 


-0.339+0.059 


8 


UGC 12545 


23 21 41.9 +27 04 15 


-20.01 


HII 


SBcd 


21.9 


8.88 


11.81 


-0.381+0.157 


9 


UGC 12546 


23 21 41.2 +27 05 14 


-20.03 


AMB 


Sbc 


21.9 


8.91 


10.64 


-0.421+0.094 



a Spectral classes and central metallicities are derived from nuclear long-slit spectral fluxes from [Barton et alj j2000fl . We use the IKewlev et al.1 
(2006b) spectral classification scheme and the Kcwlcy & Dopita (2002) [NII]/[OII] mctallicity calibration. 

Hubble Types are obtained from the NASA Extragalactic Database. 
c Separation refers to the projected separation of the pair or projected distance to the closest galaxy in the T^tuple. 
^ Gradient slope is defined in terms of ^Arf/n^^ 1 where R25 is the B-band isophote at a surface brightness of 25 mag arcsecond~ 2 . 



equal-mass disk galaxy merger models with a range of 
initial geometries and initial pericenters. We find that 
between first and second pericenter (i.e. ~ 1 Gyr), the 
metallicity gradient becomes significantly flattened and 
that the central metallicity becomes diluted by infalling 

low metallic ity gas. 

Recently, iMontuori et~aTl (|2010D used major merger 
simulations to show that the circumnuclear metallic- 
ities become diluted after the first pericenter pas- 
sage due to major gas inflows. These model predic- 
tions a r e consistent with o bservations by IKewlev et al.1 
(2006a); Rup ke et all (|2008l ) and with the simulations of 
Rupk e et al.l (|2010a|) . Montuori et al. find a second gas 
dilution peak at final coalescence in major mergers and 
significant central metallicity dilution in galaxy fly-bys, 
highlighting the importance of the gas-phase metallicity 
for tracing the recent merger or interaction history of 
galaxies. 

To date, metallicity gradients have been investigated 
for relatively small numbers of star clusters and H II re- 
gions in a small number of interacting galaxies. In this 
paper, we present the first systematic investigation of 
metallicity gradients in interacting galaxies. We present 
Keck LRIS spectra for H II regions in 8 galaxies in 
close pairs selected on the basis of their position on the 
luminosity-metallicity relation. We show that all close 
pair galaxies have significantly flatter metallicity gradi- 
ents than metallicity gradients measured in isolated spi- 
ral galaxies. We suggest that flat metallicity gradients 
provide a 'smoking gun' for recent major gas flows in in- 
teracting galaxies. Throughout this paper, we adopt the 
flat A-dominated cosmology as me asured by the WMA P 
experiment (h = 0.72, D, m = 0.29; iSpergel et al.ll2003l) . 

2. SAMPLE SELECTION AND OBSERVATIONS 

We selected 5 luminous (Mb < —20) sets of close pairs 
(separ ation between 15 - 25 kpc) from the iBarton et al.l 
(2000) galaxy pairs sample. Our pairs were selected 
to (a) span a range of positions on the luminosity- 
metallicity relation in IKewlev et al.l (|2006al ) , (b) contain 
at least one galaxy with suitable orientation for metallic- 
ity gradient analysis, and (c) to span a range of galaxy 
properties that may influence metallicity gradients (pair 
separation, presence of bars, and Hubble type). In Ta- 
ble 1, we list our sample galaxies and their properties. 



We obtained spectra for 12 - 40 star- forming regions in 
8 of our close pair galaxies with the Keck Low- Resolution 
Imaging Spectrograph (LRIS). Slitmasks were designed 
using 1" wide slits based on Ha images of the target 
galaxies observed with the Vatican Advanced Technology 
Telescope (Barton, E. & Jansen, R., private communica- 
tion). We used the 560nm dichroic with the 900 1/mm 
red grating and the 600 1/mm grism, giving a spectral 
resolution of 3.5A FWHM between 3500 - 5500A and 
2.8A between 5500 - 7000A. 

The spectra were reduced using a custom pipeline de- 
veloped in IRAF and IDL that includes standard bias 
subtraction, flatfielding, rectification, and wavelength 
calibration. Spectra were flux calibrated using standard 
stars observed at a similar airmass and time to the slit- 
mask observations. We derived emission-line fluxes from 
each spectrum using UHSPECFIT, a spectral line-fitting 
routine developed by us. F or continuum fitt ing, UH- 
SPECFIT utilizes code f rom iMoustakasI (|2006l ) and the 
iGonzafez Delgado et al.l (|2005l ) evolutionary stellar pop- 
ulation synthesis models. 

Emis sion-line fluxes were corrected for extinction us- 
ing the iCardelli et al.l (|1989l ) extinction curve. The ma- 
jority of our H II region spectra do not have the elec- 
tron temperature-sensitive auroral [O III] A4363 line, 
which is only produced in high electron temperature, low 
metallicity regions. We derived metallicities [in units of 
log(0/H) + 12] using the strong emission lines. For the 
[N II] A6584/[0 II] A3727 ratios covered by our sam- 
ple [log([N II] /[O II])> -1.3; log(0/H) + 12> 8.4], the 
[N II] /[O II] ratio is a strong function of metallicity with 
very little (<< 0.1 dex) depe ndence on the ion i zation 
state of the gas. We use the IKewlev fc Dopital (|200l 
[N II] /[O II] calibration to derive metallicities for the 
H II regions in our sample, after removal of spectra po- 
tentially contaminated by an AGN or other no n-thermal 
sources according to the (Kewlev et al.l (|2006bD classifica- 
tion scheme. More details on the spectroscopic pipeline, 
UHSPECFIT, emission-li ne flux measurement s and clas- 
si fication are provided inlRupke et all (|2010b[ ). 

IKewlev &; Ellison! (|2008| ) show that large discrepancies 
exist among metallicities derived using different strong- 
line calibrations. These discrepancies are systematic, 
therefore relative metallicities such as gradients may be 
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Fig. 1. — The B-band luminosity-metallicity (L-Z) relation from Kcwlcv ct al. (2006a) for the close pair members (red) compared with 
isolated field galaxies from the Nearby Field Galaxy Survey (unfilled circles). All metallicities in the L-Z relation are calculated from 
spectroscopy with a central long-slit aperture that covers approximately 10% of the galaxy B-band light (~ 1 kpc). The 8 close pairs 
with metallicity gradients (Table 1) are indicated on the L-Z relation as star forming galaxies (circle outlines) and composite or ambiguous 
galaxies (square outlines). Top and right panels indicate the location and metallicity of the star- forming regions in each galaxy on Ha 
images. Metallicity colors are shown in the legend bar in units of log(0/H) + 12 and the coordinates for each panel are given in units of 
arcseconds. Unfilled white circles indicate the galaxy center determined using 2MASS images. 



reliably compared if the same calibration is used. In ad- 
dition, we have verified that we obtain the same results 
if independent metallicity calibrations based on different 
line ratios are applied such as the [O HI] /H0 / [N II] /Ha 
calibration by Pcttini fc Pagell (120041). an d th e R23 cali- 
bratio ns by iKobulnickv fc Kewlevl (|2004f ) and IMcGaughl 
(p9ll) . 

3. ANALYSIS 

In Figure[T] we show the position of our sample pairs on 
the luminosity-metallicity relation (where possible) and 
the spatial distribution of metallicities for each galaxy 
(colored circles). Galaxy ID numbers are given in Ta- 
ble [1] Spatial distributions of metallicities range from 
flat (e.g., NGC 3994, #6) to relatively steep towards the 
central regions (e.g., UGC 12546, #9). 

We measure metallicity gradients for the 8 pair galax- 
ies with at least 10 H II region metallicity measurements. 
Galactocentric radii are computed using the direct dis- 
tance to the galaxy center, where the galaxy center is 
defined using 2MASS images from a well-defined bulge. 



Our distances take inclination into account using the 
optical diameters, inclination s, and line-of-node s posi- 
tion angles from HyperLeda (jPaturel et al.l [20 03) . with 
double-checks on the position angles by visual inspection 
of HI maps, where available. Figure [2] shows the metallic- 
ity gradients (in R/R25) for the close pair galaxies in our 
sample. The H II region radius is given in units of R/R25 
where R25 is the B-band isophote at a surface brightness 
of 25 mag arcsecond~ 2 . 

We directly compare the metallicity gradients mea- 
sured for our cl o se pa irs in Figure [3] for both the 
IKewlev k Dopital (l200l) [N II] /[O II] and IMcGaughl 
( 19911 ) R93 calibrations. For comparison, solid lines in- 
dicate the metallicity gradients measured in the isolated 
spiral galaxies M101, M83 and the Milky Way. For con- 
sistency, we recalculate the H II regi on metallicities in 
the is olated spiral gal axies using the IKewlev fc Dopital 
([2003) [N II] /[O II] and IMcGaughl tl991|) ca librations us- 
ing th e H II region e mission-line fluxe s fromlShaver et al.l 
(fl983T) [Milky Wavl. l Kennicutt et al.l (120031) [M101], and 
iBresolin et all (J2005] |2009j) [M83]. The metallicity gra- 



4 





NGC 3994 

0.2 0.4 0.6 0.8 



8.4 
9.1 

9. 



UGC 312 
0.2 0.4 0.6 0.8 1.0 1.2 




NGC 3995 

0.5 1.0 1.5 





0.5 1.0 
R/R 25 



0.2 0.4 0.6 0.8 1.0 
R/R 25 



Fig. 2. — Metallicity gradients derived for our close pair galax- 
ies. Solid lines indicate the least squares line of best fit to the 
metallicity gradient and the dotted lines represent the average rms 
error (0.05 dex) about the line of best fit. The metallicity gradi- 
ents are given in terms R/R25 where R is the radius and R25 is 
the radius of the B-band isophote at a surface brightness of 25 mag 
arcsecond - 2 . Similar gradients are obtained using radius in kpc. 
The same gradients are obtained withi n the errors if alter native 
metallicity calibrations are applied (sec Rupkc ct al. 2010b, for a 
detailed comparison). 

clients in these comparison galaxies ar e representative of 
isolat ed spiral galaxies (see review by IHenrv fc Worthevl 
1999). We note that we obtain the same results relative 
to the comparison galaxies if distance between the H II 
location and the galaxy center is defined in kiloparsecs. 

The metallicity gradients of all of our close pairs are 
flatter than the gradients observed in the three iso- 
lated spiral comparison galaxies, regardless of calibra- 
tion used. The mean gradient in our close pairs is 

A ar/r H ^ = — 0-25 ± 0.02 compared with a mean gra- 

dient 0f A A R/R{f } = ~° - 67 ± °' 09 for tllC isolated s P iral 

galaxies. In lRupke et al.l (|2010bl ) . we show that the mean 
metallicity gradient for a large sample of 11 isolated spi- 
ral galaxies is A A°r? = —0.57 ± 0.05. These results 
ak/k 25 

provide the first direct evidence that the metallicity gra- 
dients in galaxy pairs are systematically different from 
gradients in isolated spiral galaxies. Recent infall of less 
enriched gas from the outskirts of the close pair galaxies 

is th e likely culpri t. 

In iRupke et al.l (|2010aD . we investigate the effect of 
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Fig. 3. — Metallicity gradients for the 8 close pair galaxies in 
our sample that have metallicities measured for > 10 HII regions. 
Galaxies are labelled according to Table 1. All metallicity gradients 
are least-squares fits to HII region metallicities measured using the 
I Kewlev fc D opita ( 2002) [Nil] /[Oil] calibration (top panel) and the 
IMcGau gh (1991) R23 calibration (bottom panel). For comparison, 
we show the metallicity gradients for three isolated spiral galaxies: 
Milky Way (black), M83 (red), and M101 (blue). The 8 close pairs 
have substantially flatter metallicity gradients than the isolated 
spiral galaxies, providing a smoking gun for major merger-induced 
gas flows. 



gas inflows on metallicity gradients in theoretical merger 
simulations. We assume the simplest conditions: no on- 
going star formation, equal-mass progenitors and small 
gas mass fractions. In these theoretical models, we find a 
dramatic flattening of the initial radial metallicity gradi- 
ent between first and second pericenter. This flattening 
reflects the effects of gas redistribution over the galaxy 
disks, including both metal-poor inflow and the growth 
of tidal tails that carry metals out to large radius. 

If large merger-triggered gas flows drive the bulk of 
the inflowing gas into the central kiloparsec, triggering a 
burst of star formation, a recent nuclear burst of star for- 
matio n may be associated wi t h a fl atter metallicity gra- 
dient. iBarton Gillespie et al.l (|2003T ) applied stellar pop- 
ulation synthesis models to the B-R colors and Balmer 
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equivalent widths of the galaxy pairs sample to derive the 
fraction of R-band light originating from a recent burst 
of star formation, term ed the central burst strength. 
In iKewlev etafl (|2006al ). we showed that central burst 
strength correlates strongly with central metallicity in 
close pairs. T his result is consi s tent w ith the theoretical 
simulations of Montuo ri et al.l ([2010), which show that 
pairs experiencing the strongest bursts of star formation 
also have the strongest circumnuclear dilution. If most 
of the infalling gas is deposited in the central regions, 
we might expect a correlation between the metallicity 
gradient slope and central burst strength in our galaxy 
pairs. However, we find no statistically significant cor- 
relation between the metallicity gradient and the central 
burst strength (Spearman Rank coefficient of -0.14 with 
a two-sided probability of obtaining this value by chance 
of 74%). Similarly, we find no statistically significant cor- 
relation between the metallicity gradient slope and the 
presence of blue bulges or pair projected separation. We 
note that our gradient sample is limited by the small 
sample size and that a larger sample is required for more 
robust conclusions on the relationship between the metal- 
licity gradient slope and galaxy properties. If the lack of 
correlation between central burst strength and metallic- 
ity gradient slope holds for a larger sample, these results 
would suggest that (a) the timescale for a central star- 
burst differs from the timescale over which the metallic- 
ity gradient remains disrupted, and/or (b) the metallicity 
gradient slope is strongly influenced by radial gas redis- 
tribution along the spiral arms and/or metal outflows, 
rather than a simple inflow into the central kiloparsec. 

The presence of bars is correlated with efficient gas 
flows to ward the nucleus in late-type spiral g alax- 
ies (e.g., iSakamoto et~aTl 119991 : iSheth et all 120051 ) and 
bars may be created in t he early stages of a merger 
([Barnes fc Hcrnquist 199ll). Two galaxies in our sample 
have bars (UGC 312 and UGC 12545). While UGC 312 
has a low central metallicity and relatively flat gradient 

-0.133), UGC 12545 has a steeper 



gradient slope ( A a°r/r^ 



-0.381) that is within the 



full range of gradient slopes observed in our pairs sam- 
ple. It is not clear from these two galaxies whether bars 



allow more efficient transport of gas towards the central 
regions during galaxy mergers. 

We find no correlation between the distance from the 
nearby field galaxy LZ relation and the close pair metal- 
licity gradient slope. Galaxy pairs that lie close to 
the field galaxy LZ relation in Figure 1 may have had 
high central metallicitie s [log(0/H) + 12> 9.1 in the 
IKewlev fc Dopital (120021) [N Hl/[0 II] scale] prior to the 
interaction. In iRupke et ail ([2010aD . we investigated the 
theoretical relationship between the metallicity gradient 
and central metallicities. We showed that metallicity gra- 
dients become disrupted by the interaction, followed by 
a dilution in the central metallicity. The observational 
results presented here support this scenario. 

4. CONCLUSIONS 

We present the first systematic study of metallicity gra- 
dients in close galaxy pairs. We find that all of our galaxy 
pairs have flatter metallicity gradients than gradients ob- 
served in typical isolated spiral galaxies. These results 
show that the low central metallicities seen previously in 
galaxy pairs and interacting galaxies signal a relatively 
flat metallicity gradient. We conclude that large gas in- 
flows are responsible for the disruption and flattening of 
metallicity gradients in close pair galaxies. Evidence for 
gas flows has been difficult to study in the past, requiring 
integral field or neutral gas velocity maps. Our results 
indicate that metallicity gradients provide a unique and 
feasible method for detecting recent strong tidal gas flows 
in galaxy interactions and mergers. Our future work 
includes comparisons between metallicity gradients and 
theoretical merger simulations and an extension of our 
pairs sample to later stage mergers. 



We thank the referee for comments that improved this 
paper. We thank Rolf Jansen for providing Ha data used 
for slitmasks and Li Hsin Chien for preliminary work on 
some of this data. The authors thank the indigenous 
Hawaiian community for use of their sacred mountain. 
This research has made use of NASA's Astrophysics Data 
System Bibliographic Services and the NASA/IPAC Ex- 
tragalactic Database (NED). 



REFERENCES 



Barnes, J. E., & Hernquist, L. 1996, ApJ, 471, 115 
Barnes, J. E., & Hcrnquist, L. E. 1991, ApJ, 370, L65 
Barton, E. J., Geller, M. J., k. Kenyon, S. J. 2000, ApJ, 530, 660 
Barton Gillespie, E., Geller, M. J., & Kenyon, S. J. 2003, ApJ, 
582, 668 

Bastian, N., Trancho, G, Konstantopoulos, I. S., & Miller, B. W. 

2009, ApJ, 701, 607 
Bresolin, F., Ryan- Weber, E., Kennicutt, R. C, & Goddard, Q. 

2009, ApJ, 695, 580 
Bresolin, F., Schaerer, D., Gonzalez Delgado, R. M., & Stasinska, 

G. 2005, A&A, 441, 981 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 

245 

Chien, L., Barnes, J. E., Kewley, L. J., & Chambers, K. C. 2007, 
ApJ, 660, L105 

Cooper, M. C, Tremonti, C. A., Newman, J. A., & Zabludoff, 

A. I. 2008, MNRAS, 390, 245 
Cullen, H., Alexander, P., Green, D. A., & Sheth, K. 2007, 

MNRAS, 376, 98 
Ellison, S. L., Patton, D. R., Simard, L., & McConnachie, A. W. 

2008, AJ, 135, 1877 



Ellison, S. L., Simard, L., Cowan, N. B., Baldry, I. K., Patton, 
D. R., &; McConnachie, A. W. 2009, MNRAS, 396, 1257 

Emonts, B. H. C, Morganti, R., Tadhunter, C. N., Holt, J., 
Oosterloo, T. A., van der Hulst, J. M., & Wills, K. A. 2006, 
A&A, 454, 125 

Georgakakis, A., Forbes, D. A., & Norris, R. P. 2000, MNRAS, 
318, 124 

Gonzalez Delgado, R. M., Cervino, M., Martins, L. P., Leitherer, 

C, & Hauschildt, P. H. 2005, MNRAS, 357, 945 
Henry, R. B. C, & Worthey, G. 1999, PASP, 111, 919 
Hibbard, J. E., & van Gorkom, J. H. 1996, AJ, 111, 655 
Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., 

Robertson, B., & Springel, V. 2006, ApJS, 163, 1 
Horellou, C, & Koribalski, B. 2007, A&A, 464, 155 
Iono, D., Yun, M. S., & Ho, P. T. P. 2005, ApJS, 158, 1 
Kennicutt, R. C, Bresolin, F., & Garnett, D. R. 2003, ApJ, 591, 
801 

Kewley, L. J., & Dopita, M. A. 2002, ApJS, 142, 35 
Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183 
Kewley, L. J., Geller, M. J., & Barton, E. J. 2006a, AJ, 131, 2004 
Kewley, L. J., Groves, B., Kauffmann, G, & Heckman, T. 2006b, 
MNRAS, 372, 961 



6 



Kobulnicky, H. A., & Kewley, L. J. 2004, ApJ, 617, 240 
Lee, J. C, Salzer, J. J., & Melbourne, J. 2004, ApJ, 616, 752 
Martin, C. L. 2006, ApJ, 647, 222 
McGaugh, S. S. 1991, ApJ, 380, 140 

Michel-Dansac, L., Lambas, D. G., Alonso, M. S., & Tissera, P. 

2008, MNRAS, 386, L82 
Montuori, M., Di Matteo, P., Lehnert, M. D., Combes, F., & 

Semelin, B. 2010, ArXiv e-prints, 1003.1374 
Moustakas, J. 2006, PhD thesis, The University of Arizona, 

United States - Arizona 
Narayanan, D. et al. 2008, ApJS, 176, 331 

Paturel, G., Petit, C, Prugniel, P., Theureau, G., Rousseau, J., 
Brouty, M., Dubois, P., & Cambresy, L. 2003, A&A, 412, 45 
Peeples, M. S., Pogge, R. W., & Stanek, K. Z. 2009, ApJ, 695, 259 
Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 
Rampazzo, R., Plana, H., Amram, P., Bagarotto, S., Boulcstcix, 

J., & Rosado, M. 2005, MNRAS, 356, 1177 
Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005, ApJS, 160, 87 



Rupke, D. S. N., Kewley, L. J., Chien, L.-H. 2010b, ApJ, 00, 
submitted 

— . 2010a, ApJ, 710, L156 
Rupke, D. S. N., Veilleux, S., & Baker, A. J. 2008, ApJ, 674, 172 
Sakamoto, K., Okumura, S. K., Ishizuki, S., & Scoville, N. Z. 

1999, ApJS, 124, 403 
Shaver, P. A., McGee, R. X., Newton, L. M., Danks, A. C, & 

Pottasch, S. R. 1983, MNRAS, 204, 53 
Sheth, K., Vogel, S. N., Regan, M. W., Thornley, M. D., & 

Teuben, P. J. 2005, ApJ, 632, 217 
Spergel, D. N., et al. 2003, ApJS, 148, 175 

Springel, V., Di Matteo, T., & Hernquist, L. 2005, ApJ, 620, L79 
Sol Alonso, M. S., Michel-Dansac, L., & Lambas, D. G. 2010, 
A&A, 514, 57 

Trancho, G., Bastian, N., Miller, B. W., & Schweizer, F. 2007a, 
ApJ, 664, 284 

Woods, D. F., Geller, M. J., & Barton, E. J. 2006, AJ, 132, 197 



